Self-oscillations in a superconducting stripline resonator integrated with a DC-SQUID 
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We study self-sustained oscillations (SO) in a Nb superconducting stripline resonators (SSR) 
integrated with a DC superconducting quantum interface devices (SQUID). We find that both the 
power threshold where these oscillations start and the oscillations frequency are periodic in the 
applied magnetic flux threading the SQUID loop. A theoretical model which attributes the SO 
to a thermal instability in the DC-SQUID yields a good agreement with the experimental results. 
This flux dependant nonlinearity may be used for quantum state reading of a qubit-SSR integrated 
device. 



We study thermal instability in superconducting 
stripline resonators (SSR's) working at gigahertz fre- 
quencies. We have recently demonstrated how thermal 
instability can create extremely strong nonlinearity in 
such resonators^. This nonlinearity is manifested by 
self-sustained oscillations (SO) at megahertz frequencies, 
strong intermodulation gain, stochastic resonance, sensi- 
tive radiation detection, and more^. In the present work 
we have integrated a SSR with a DC superconducting 
quantum interface devices (SQUID). Similar configura- 
tions have been recently studied by other groups and 
it was shown that such devices can be used as readout 
and coupling elements for qubits^i^. In addition, it was 
shown that further improvement in the sensitivity of the 
readout process of the qubit is possible by biasing the 
resonator to a state of nonlinear responsitivit y 5 ' 6 ' 7 ' 8 ' 9 '. 

Our experiments are performed using the setup de- 
picted in Fig. QJa). We study the response of the in- 
tegrated SSR to a monochromatic injected pump tone 
that drives one of the resonance modes, and measure 
the reflected power spectrum by a spectrum analyzer. 
We find that there is a certain range in the plane of 
the pump-frequency pump-power parameters, in which 
intrinsic SO occur in the resonator. These oscillations 
are manifested by the appearance of sidebands in the re- 
flected power spectrum. In addition, we apply bias mag- 
netic flux through the DC-SQUID and find that both 
the threshold where these oscillations start and their fre- 
quency are periodic in the applied magnetic flux, hav- 
ing a periodicity of one flux quantum. We extend our 
theoretical model^, which has originally considered SO 
and thermal instability in a SSR integrated with a single 
micro-bridge, to include this flux dependency and find a 
good agreement with the experimental results. 

A simplified circuit layout of our device is illustrated 
in Fig. DJb). We fabricate our devices on a Silicon wafer, 
covered by a thin layer of Silicon Nitride. Each device 
is made of a thin layer (< 100 nm) of Niobium and com- 
posed of a stripline resonator having a DC-SQUID Jljc)) 
monolithically embedded into its structure. The res- 
onator length is I = 18 mm, and its first resonance mode 
is found at /i = cji/2tt = 3.006 GHz. The DC-SQUID 
has two nano-bridges dd)), one in each of its two arms. 
Their size is typically 100 nm 2 and therefore each nano- 
bridge functions as a weak-link that approximately can 




FIG. 1: (a) Measurement setup, (b) Schematic layout of our 
device, (c) SEM image of the DC-SQUID, (d) SEM image of 
the two nano-bridge-based Josephson junctions. 



be regarded as a regular Josephson junctio n 10 ! 11 . A feed- 
line, weakly coupled to the resonator, is employed for de- 
livering the input and output signals. An on-chip filtered 
DC bias line passing near the DC-SQUID is used to apply 
magnetic flux through the SQUID. Some measurements 
are carried out while the device is fully immersed in liq- 
uid Helium, while others in a dilution refrigerator where 
the device is in vacuum. Further design considerations 
and fabrication details can be found elsewhere^. 

Figure 2 exhibits SO in the power reflected off the 
resonatorii^. In this measurement we inject into the 
resonator an input pump tone having a monochromatic 
frequency oj v = oji, and measure the reflected power 
spectrum around u>\ while varying the input pump 
power Pp. At relatively low (P p < — 66dBm) and high 
(P p > — 48dBm) pump powers the response of the res- 
onator is linear, namely, the reflected power spectrum 
contains a single spectral component at the frequency of 



2 




FIG. 2: Typical experimental results of the SO phenomenon. 
The color-map shows the reflected power P rc fl as a function of 
the input pump power P p and the measured frequency /sa, 
centered on the pump frequency, which coincides with the 
first resonance frequency /i (AfsA ~ fsA — Note that 
the data is truncated at P re fl = —35 dBm for clarity. The 
inset shows a cut of that measurement obtained with P p — 
-61dBm. 



the stimulating pump tone. In between these two power 
thresholds the resonator self oscillates and regular modu- 
lation of the pump tone occurs. As a result, the reflected 
spectrum contains several orders of modulation products 
realized by rather strong and sharp sidebands (see inset 
of Fig. [2]) that extend to both sides of the pump tone 
frequency. The SO frequency, defined as the frequency 
difference between the pump frequency and the primary 
sideband, ranges between few to tens of megahertz and 
increases with the pump power. 

The dependence of the SO on the applied mag- 
netic flux is shown in Fig. [3J In this measurement 
we inject a pump tone having a stationary frequency 
(u) p ~ oji) and power (P p = — 51.5dBm), and measure 
the reflected power spectrum around lu± while vary- 
ing the DC bias current. As shown, the SO fre- 
quency is periodically changed by the applied magnetic 
flux. The periodicity is one flux quantum and the rel- 
ative change is typically about 20%. Figures HJa) and 
(c) show two measurements of flux dependent SO ob- 
tained with pump power above and equal the thresh- 
old power (-P p ,( a ) = — 51.5dBm, P p j c ) — — 51.6dBm), re- 
spectively. The later measurement demonstrates how 
magnetic flux can switch the resonator from a steady 
state response to a limit-cycle state where it experiences 
SO. 

To account for our results we model our device as a 
transmission line resonator interrupted by a DC-SQUID. 
The impedance of the DC-SQUID is composed of a con- 
stant inductor in series with a flux-dependent inductor 
shunted by a parallel resistor —. The flux dependence of 
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FIG. 3: Typical experimental results of flux dependent SO. 
The color-map shows the reflected power P re fl as a function 
of the bias DC current Jfl ux and the measured frequency /sa 
centered on the pump frequency / p = /i (A/sa = /sa — /i). 



the DC-SQUID impedance gives rise to periodic depen- 
dence of the SSR's damping rates on the applied mag- 
netic field (the change in the SSR resonance frequency is 
relatively small). 

The dynamics of our system can be captured thus by 
two coupled equations of motion^. The first equation 
describes the dynamics of a resonator driven by feed-line 
carrying an incident coherent tone b m = & n e _lWpt , where 
6g n is constant complex amplitude and ui p = 2nf p is the 
driving angular frequency. The mode amplitude inside 
the resonator can be written as _Be _lWpt , where B (t) is 
complex amplitude, which is assumed to vary slowly on 
a time scale of l/w P - In this approximation, assuming a 
noiseless system, the equation of motion of B readsi^ 

dB 

— = [i(u; p -u; )- 1 ]B~iy/2T 1 b i », (1) 

where ujq is the angular resonance frequency and 7 (T) = 
7i + 72 (T, &ex), where 71 is the coupling coefficient be- 
tween the resonator and the feed- line and 72 (T, <& ex ) is 
the temperature and flux dependent damping rate of the 
mode. 

The heat-balance equation for the temperature T of 
the nano-bridges composing the DC-SQUID is given byA2, 

C^- =2^72 \B\ 2 -H{T- T ), (2) 
dt 

where C is the thermal heat capacity, H is the heat trans- 
fer coefficient, and To is the temperature of the coolant. 

Coupling between Eqs. ([T]) and ^ originates by the 
dependence of the damping rate 72 (T, <& ex ) of the driven 
mode on the impedance of the DC-SQUID 15 , which in 
turn depends on the temperature of its nano-bridges and 
on the applied magnetic flux. We assume a simple case, 
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FIG. 4: Typical experimental and numerical results of flux 
dependent SO. Each panel shows a color-map of the reflected 
power P re fl as a function of the bias DC current 7fl ux and 
the measured frequency /sa centered on the pump frequency 
f p — fi (A/sa = fsA — fi)- Panels (a) and (c) show experi- 
mental results obtained at pump power above and equals the 
low threshold power, respectively. Panels (b) and (d) show 
the corresponding theoretical results obtained by numerically 
integrating the equations of motion. 



where the dependence on the temperature T is described 
by a step function that occurs at the critical temperature 
T c of the superconductor. We further assume that only 
when the nano-bridges are in a superconducting phase 
the damping rate depends on the external flux. Namely, 



72 takes a flux dependent value 72s ($ e i) when the nano- 
bridges are in a superconducting phase (T < T c ) and a 
flux-independent value 72 n when they are is a normal- 
conducting phase (T > T c ). 

The results of a numerical integration of the equations 
of motion are shown in Fig. |4fb) an d (d). Our model 
qualitatively reproduces the same flux dependency of the 
SM frequency on the magnetic flux. The parameters 
used for the numerical simulation were obtained as fol- 
lows. The thermal heat capacity C = 4.5 nJ cm~ 2 K -1 
and the heat transfer coefficient H = 24.5 mWcm~ 2 K _1 
were calculated analytically according to Refs i 16 ' 17 . The 
values of the resonance frequency and the various damp- 
ing rates, 7i s = 2.4 MHz, 72s £ [2.93,3.32] MHz, 72n = 
13 MHz, and uiq — 2ir ■ 3.006 GHz were extracted from 
Sll reflection coefficient measurements according toi£. 

In conclusion, we report on periodic flux-dependency 
of SO in a SSR integrated with a DC-SQUID. The flux 
significantly modulates the oscillation frequency and can 
be used to turn them on and off if the device is driven 
near the power threshold. A theoretical model which 
attributes this behavior to thermal instability in the DC- 
SQUID exhibits a good quantitative agreement with the 
experimental results. 
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